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ANALYSIS OF LOSS-OF-COOLANT ACCIDENT FOR A FAST-SPECTRUM LITHIUM- 
COOLED NUCLEAR REACTOR FOR SPACE-POWER APPLICATIONS 
by George E. Turney, Edward J. Petrik, and Arthur W. Kieffer 

Lewis Research Center 

SUMMARY 

A two-dimensional, transient, heat-transfer analysis was made to determine the 
temperature response in the core of a conceptual space-power nuclear reactor following 
a total loss of reactor coolant. With loss of coolant from the reactor, the controlling 
mode of heat transfer is thermal radiation. 

In one of the schemes considered for removing decay heat from the core, it was 
assumed that the 4n shield which surrounds the core acts as a constant- temperature 
sink (temperature, 700 K (1260° R)) for absorption of thermal radiation from the core. 
Results based on this scheme of decay heat removal show that the surface emissivity 
has a significant effect on the temperatures in the core. With an emissivity of 0. 20, 
the core centerline fuel pin reached the melting temperature of UN (2923 K, or 5260° R) 
at about 1900 seconds after the loss- of- coolant accident. And with emissivity values of 
0. 40 and 0. 80, the maximum fuel-pin temperatures were 2820 K (5080° R) and 2355 K 
(4240° R), respectively. Thus for emissivity values greater than about 0. 40, thermal 
radiation to only the surrounding 4 n shield is sufficient to prevent the fuel from reach- 
ing its melting point following a loss of coolant. 

In another scheme for removing decay heat, the centerline fuel pin was replaced by 
a redundant coolant channel. With this arrangement, decay heat is transferred to both 
the redundant coolant channel and the 4n radiation shield. A constant temperature of 
1222 K (2200° R) was assumed for the redundant coolant channel. Based on an emissiv- 
ity value of 0.20 for all core materials, the calculated maximum fuel temperature for 
this scheme of heat removal was 2840 K (5100° R), or about 90 K (160° R) less than the 
melting temperature of uranium nitride (UN). Hence, it appears that a single coolant 
channel at the center of the core is sufficient to prevent fuel melting, even with an emis- 
sivity value as low as 0. 20. 


INTRODUCTION 


Our nation’s future space missions will require systems capable of delivering large 
amounts of electric power. To meet these future requirements, the Lewis Research 
Center has been working on a technology program aimed at the design of an advanced, 
high-powered, nuclear Brayton space powerplant. The heat source being considered 
for this advanced space powerplant is a compact, fast- spectrum, nuclear reactor. The 
basic design requirements for the reactor are that it produce 2. 17 megawatts thermal 
power and that it operate continuously for 50 000 hours. In this conceptual space power- 
plant, the heat generated by the reactor is removed by liquid lithium which circulates 
continuously through a closed primary flow loop. The primary loop of this system is 
thermally coupled to one or more complete inert- gas Brayton power conversion loops. 
Neutronic design information for this fast- spectrum reactor is presented in reference 1; 
and the operating characteristics of the primary flow loop of this system are described 
in reference 2. 

An important area which must be considered in the overall design of this system is 
the analysis of potential failures and/or system malfunctions. A preliminary analysis 
of several possible accidents in this system is reported by Davison (ref. 3). Of the 
many types of malfunctions and/or failures possible, one of the most serious is that 
which results from a complete loss of coolant in the reactor primary flow loop. With 
loss of coolant from the core, all heat interchange between the fuel, interior core struc- 
ture, reflector, pressure vessel, and reactor shield is by radiation. For a failure of 
this type, we would like to know the temperature response of the fuel pins inside the 
core and also the conditions under which melting of fuel in the core is possible. 

Two separate schemes for removing heat from the core were considered. In one 
scheme, we assumed that the radiation shield which surrounds the core acts as a 
constant- temperature sink for the absorption of heat transferred from the core. In the 
other scheme, the centerline fuel pin in the core was replaced by a redundant coolant 
channel. With this arrangement, there are two sinks for absorption of reactor heat: 
one at the center of the cylindrically shaped core, and the other outside the core (i. e. , 
the surrounding radiation shield). 

In order to investigate the loss- of- coolant accident for these two schemes of heat 
removal, a two-dimensional, transient, heat- transfer analysis was made. From the 
analysis, we determined the temperature response of the fuel pins inside the core and 
the transient temperature distribution throughout the core. The analysis was made by 
using the Chrysler Improved Numerical Differencing Analyzer for 3rd Generation Com- 
puters (CINDA-3G) digital computer program (ref. 4). The results of this study and a 
description of the thermal network model used to represent the core are presented in 
this report. 
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DESCRIPTION OF REACTOR 


The reference design for this fast- spectrum reactor has 247 cylindrical fuel pins, 
each with an outside diameter of 1. 905 centimeters (0. 750 in. ) and a length of 37. 6 cen- 
timeters (14. 8 in. ). Figure 1 is a schematic of the reactor core design. 

The main body of the core has 181 fuel pins arranged symmetrically in a six-pointed 
star shape. Another 66 fuel pins are in the six control drums just outside the main body 
of the core. Figure 2 shows a cross-sectional view of the fast- spectrum reactor core. 

The fuel pins in the core contain cylindrically shaped uranium nitride (UN) pellets 
which are covered with a tungsten liner and encased in a T-lll (Ta-8W-2Hf) cladding. 
The interior of each fuel pellet has a central void which provides space for fission gas 
containment and fuel swelling. Figure 3 shows the arrangement of UN pellets in the 
present fuel- pin design. 

The fuel pins in the core are cooled by lithium which flows through annular passages 
formed by the outside surfaces of the fuel pins and the inside surfaces of concentric 
T-lll tubes which surround the fuel pins. The surrounding T-lll tubes form a honey- 
comb structure; the tubes in this structure have an inside diameter of 2. 11 centimeters 
(0. 83 in. ) and a wall thickness of 0. 025 centimeter (0. 010 in. ). 

The reactor thermal power is regulated by six control drums which, when rotated, 
move fuel (or a T-lll neutron absorber) near or away from the main body of the core. 
Each control drum has 11 fuel pins, a reflector of TZM (Mo-0. 5Ti-0. 08Zr), and a T-lll 
neutron absorber. Between the control drums near the periphery of the core is a fixed 
reflector of TZM. The pressure vessel surrounding the core is T-lll with a wall thick- 
ness of 0. 625 centimeter (0. 250 in. ). 

In the conceptual system design, the entire pressure vessel is surrounded with a 
radiation shield (477 type) consisting of alternate layers of tungsten and lithium-6 hydride. 

As indicated in figure 2, the fuel pins are grouped into three zones. Fuel enrich- 
ment for each zone is the same, that is, 93.2 atomic percent of the uranium is the 235- 

nor 

isotope (U^ ja ). However, the diameter of the central void in the fuel pellets (and con- 
sequently the volume fraction of UN) is different for each zone. Figure 4 shows a cross 
section of a fuel pin, along with dimensions of the pins in each of the three fuel zones. 

The design-point operating conditions of the primary flow loop which are relevant 
to this study are listed in table I. 
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Figure 2 . - Cross section of fast-spectrum, space-power reactor core. 



Figure 3. - Fuel-pin design for fast-spectrum, space-power reactor. 
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r Clad (tungsten) 



Fuel zone 

Void diameter, 

d 

Number of pins 
per zone 

cm 

in. 

I 

0. 762 

0.300 

73 

II 

.678 

.267 

90 

III 

.480 

.189 

84 


Figure 4. - Cross section of fuel pin and fuel-pin 
dimensions for the three fuel zones. 


TABLE I. - DESIGN- POINT OPERATING CONDITIONS 


Reactor thermal power, MW 2. 17 

Lithium temperature at reactor inlet, K (°R) 1167(2100) 

Lithium temperature at reactor outlet, K (°R) 1222(2200) 

Lithium flow rate, kg/sec (lbm /sec) 9.39(20.7) 
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! ANALYSIS 

| General Considerations 

If a break occurs in either the reactor pressure vessel or one of the primary- loop 
( flow lines, the entire inventory of lithium coolant would be lost from the primary loop. 

Based on the design lithium flow rate (table I) and the inventory of lithium in the core, 

I such a failure could result in the core being completely voided of lithium in less than 
| 0. 5 second. 

In all likelihood a complete loss of coolant from the primary flow loop will render 
the system permanently inoperable. That is, it is unlikely that a failure of this type 
could be corrected in space and the system returned to normal operation. Therefore, 
the primary concern associated with the loss- of- coolant accident in this space- power 
reactor is to prevent the fuel from melting and possibly recombining to form a super- 
critical assembly. 

In the event of a loss of reactor coolant, it is desirable that the reactor undergo an 
immediate (scram) shutdown. We have assumed, in the analytical treatment of this 
problem, that loss of coolant and shutdown of the reactor occur simultaneously. Hence, 
all heat generated in the core following a loss of coolant is decay heat - resulting from 
fission product decay, fissions by delayed neutrons, and radioactive emission due to 
parasitic absorption of neutrons. 

The magnitude of the decay heat generation rate in the core following shutdown de- 
pends primarily on (1) the reactor's initial operating power level and (2) the period of 
time during which the reactor was in operation before the shutdown. To prevent melting 
of fuel in the core, a significant fraction of the decay heat must be transferred by ther- 
mal radiation from the core. 

As stated in the INTRODUCTION, we investigated two separate schemes for remov- 
ing decay power from the core following a total loss of reactor coolant. In one scheme 
we assumed that the 4 jt radiation shield which surrounds the core acts as a single, 
constant- temperature sink for heat transfer from the core. With this single thermal 
sink, heat is transferred radially outward from the interior of the core to the surround- 
ing radiation shield. In the other scheme for removing decay power, the centerline fuel 
pin in the core was replaced by a redundant, constant- temperature coolant channel. With 
this arrangement, there are two thermal sinks for absorption of reactor decay heat. 

And decay heat generated in the fuel is transferred to both the redundant centerline cool- 
ant channel and the 4rr radiation shield which surrounds the core. 
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Analytical Model of Core 


The heat interchange process which takes place in the core following a loss of cool- 
ant was considered to be two-dimensional. Figure 5 is a sketch of the model used to 
represent the core for the transient heat-transfer analysis. It is a 30° sector of the core 
cross section with unit axial depth. The sides of the 30° sector model are adiabatic 
boundaries, and the control drums are assumed to be in their full- out (shutdown) posi- 
tion. 

As indicated in figure 5, the 30° sector model is divided into a number of finite vol- 
ume elements, each of which is represented by a nodal point. The boundaries of these 
volume elements are indicated by dashed lines. A total of 230 nodal points are used to 
describe the 30° sector model of the core. The breakdown of nodal points in this model 
is as follows: In the main body of the core, 96 nodes represent the fuel pins, and 
another 96 nodes represent the T-lll honeycomb tubes. In the control drum half, seven 
nodes describe the fuel pins, and 15 additional nodes describe the TZM and T-lll mater- 



Figure 5. - Model for two-dimensional heat-transfer analysis of core following 
loss of coolant. 
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ial volumes. The TZM fixed reflector region has 10 nodes, and the T-lll pressure 
shell has five nodes. The radiation shield which surrounds the core is represented by 
a single node. 

From figure 5, we see that each fuel pin in the core main body is divided by radial 
lines into a number of truncated, pie- shaped volume elements. And each of these vol- 
ume elements, consisting of UN fuel, tungsten cladding, and T-lll cladding, is repre- 
sented by a single nodal point. The individual nodal points associated with each fuel pin 
and each surrounding honeycomb tube in the main body of the core are designated by a 
number and letter. The identification scheme is illustrated in figure 5 for fuel pin 15. 
Thus, node e of fuel pin 15 represents the volume element of that fuel pin which extends 
from 9 = it to 9 = 4n/Z. Similarly, node k represents the honeycomb tube volume ele- 
ment which extends from 6 = n to 9 = 4n/3. 

For this analysis, we assumed the heat source for the fuel pins to be distributed 
uniformly throughout the UN fuel, the tungsten cladding, and the T-lll cladding. (This 
assumption was necessary since only one radial node was used for each fuel- pin volume 
element. ) Based on this assumption, the nodal points in the fuel pins were positioned 
to coincide with the location of the mean radial temperature in the fuel- pin volume ele- 
ments. The resulting radial locations of nodal points in the fuel pins of fuel zones I, II, 


TABLE n. - RADIAL LOCATION OF 
NODES IN FUEL PINS 


Fuel zone 

Nodal location (radial distance 
from fuel- pin centerline) 

cm 

in. 

I 

0.673 

0.265 

II 

.663 

.261 

III 

.640 

.252 


and m are given in table n. The nodal points in all other volume elements were 
assumed to be at the geometric center (i.e. , the centroid of area) of the elements. 

The CINDA-3G thermal analyzer code (ref. 4) uses an electrical circuit model to 
represent the thermal network of the physical problem. And the equations which des- 
cribe the thermal network are analogous to those which describe an equivalent electrical 
circuit. The thermal capacitances assigned to the nodes and the thermal resistances 
between nodes are represented by so-called ’'lumped- parameters. " These lumped- 
parameters of resistance and capacitance form a resistance- capacitance thermal net- 
work. And the set of differential equations which describe the behavior of the thermal 
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network are solved simultaneously by an iterative technique. For this study we used the 
CINDA-3G execution subroutine called CNFRWD. 

Analytical Assumptions and Input Data 

The major assumptions used in the analysis of the loss- of- coolant accident are as 
follows : 

(1) The initial temperature of the core (i.e. , the temperature of fuel, cladding, 
interior core structure, reflector, and pressure vessel) is 1222 K (2200° R). This is 
the approximate average design operating temperature of the core. 

(2) The temperature of the radiation shield which surrounds the core remains con- 
stant at 700 K (1260° R) following the loss of coolant. (The radiation shield acts as a 
sink for absorption of thermal radiation leaving the core. ) This temperature is the 
approximate average design operating value for the radiation shield. 

(3) The loss of coolant from the core is instantaneous and complete. Consequently, 
the transient phenomena associated with the loss of coolant process was not considered. 

(4) All heat interchange between material surfaces in the core is by thermal radia- 
tion. 

(5) The reactor is operated at design power for 1 year prior to the loss- of- coolant 
accident. And the reactor decay power used in the analysis was based on this design 
power operating period. 

(6) For those calculations made with a redundant coolant channel at the center of the 
core, we assumed that this channel acted as an additional sink, and its temperature re- 
mains constant at 1222 K (2200° R) following the loss-of-coolant accident. 

In order to solve the thermal network equations, we must specify the initial conditions, 
the boundary conditions, and the heating rate in each fuel node. The initial conditions 
(i. e. , the starting temperatures of all nodal points) are given by assumption (1). And 
the boundary conditions are given by assumptions (2) and (6). 

The heating rate in each fuel node decays with time after the loss of coolant. As 
stated in the section General Considerations, we have assumed in the analysis of this 
problem that the reactor is shut down concurrently with the loss of coolant and that all 
heat generated in the core following the loss of coolant is decay heat. 

Figure 6 shows the ratio of reactor decay power to design power against time after 
reactor shutdown. The data in this figure are based on the assumption that the reactor 
is operated at design power for a period of 1 year before shutdown. (As stated by 
assumption (5) in this section, the decay power used in this analysis was based on a de- 
sign power operating period of 1 year prior to shutdown. ) 

The fuel-pin power factors (i.e. , the local- to- average power ratios among the fuel 
pins) following 1 year of continuous design power operation are shown in figure 7. From 
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1_L 

Time after reactor shutdown, sec 

Figure 6. - Ratio of reactor decay power to design power against time after shutdown. Reactor operating 
time at design power before shutdown, 1 year. 



Figure 7. - Fuel-pin power factors ( local -to-average 
power ratios) following 1 year of continuous 
design power operation. 








the fuel- pin power factors (fig. 7) and the decay power data (fig. 6), the time- dependent 
heating rate for each fuel pin was determined as follows: 


'Fuel- pin heating! 
^ rate at time 9 J 


' Reactor design power \ 

^Number of fuel pins in reactory 

( Ratio of reactor decay! 

power to design power I X (Fuel- pin power factor) 
at time 9 J 


( 1 ) 


Similarity, the time- dependent heating rate in each fuel-pin node was computed as the 
product of the fuel-pin heating rate (eq. (1)) and the volume fraction of the fuel pin 
associated with the fuel- pin node. 

The thermophysical materials properties used in this study were assumed to be 
constant and independent of temperature. The values used for the specific heats and 
thermal conductivities are listed in table HI. The data in table HI were taken from the 
references indicated and were evaluated at a temperature of 1922 K (3460° R). 

Two kinds of nodes, diffusion and boundary, were used in the thermal network 
representation of the 30° sector model of the core. Each constant- temperature heat 
sink was represented by a single boundary node. All other nodes in the 30° sector 
model were diffusion nodes. Each diffusion node has a thermal capacitance associated 
with it, which permits storage of thermal energy. The thermal capacitance of a diffu- 
sion node is defined as the product of mass associated with the node and the specific 
heat (i.e., pVCp). (All symbols are defined in the appendix. ) 

Two types of thermal conductances were used in the network equations, namely 


TABLE in. - THERMOPHYSICAL PROPERTIES 


Material 

Specific heat 

Thermal conductivity 

kW-sec 

Btu 

kW 

Btu 

(kg)(K) 

(lbm)(°R) 

(m)(K) 

(sec)(ft)(°R) 

T-lll a 

17.1 

0.041 

0.0606 

0. 0097 

UN b 

27.7 

.066 

.0303 

.0049 

Tungsterf 

16.8 

.040 

.1038 

.0167 

TZM d 

29.3 

.070 | 

.0900 

.0144 


a Ref. 5. c Ref. 7. 

b Ref. 6. d Ref. 8. 
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conduction and radiation. The conductance between conductor nodes G^, is defined as 
the inverse of the thermal resistance between nodes; that is, 


G 


C 


kA 

l 


( 2 ) 


Radiation conductance does not have an equivalent inverse- resistance connotation. 
Rather, the radiation conductance between any two nodes designated as a and /3 in 
the 30° sector model is denoted by Gj^ and is defined as 


G R, a—/3 ~ aA s, a I 


+ 

A 


La 


s, a (_ 1 

S, 13 \ e t3 



Of— /3 


(3) 


The thermal capacitance, conductor conductance, and radiation conductance repre- 
sent the lumped parameters which are input to the CINDA-3G program. 


Analytical Procedure 

The equations which describe the thermal response of all nodes in the 30° sector 
model are formulated from the basic energy equation. That is, for any node designated 
a in the 30° sector model we can write 


( Rate of change in\ /Rate of net heat interchange by\ 
internal energy 1 = I conduction between node a I 
, of node a f \ and adjacent conduction nodes / 


/Rate of net heat interchange by 
+ | radiation between node a 
and adjacent radiation nodes 



Rate of heat generation in\ 
the volume represented I 
by node a J 


(4) 


In the paragraphs which follow, we describe how equation (4) is applied to determine 
the temperature response of both a typical fuel- pin node and a typical node in the T-lll 
honeycomb structure. 
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Consider, for example, a cluster of three fuel pins and honeycomb tubes with nodal 
designations, as shown in figure 8. The thermal resistance network for the fuel-pin 
node e is depicted in figure 9. With reference to figure 9, the energy balance for node 
e can be written as follows: 


AT 

p e v e C p,e"^' 


= [ G C,WVV +G C,d^< T d- T e>] 


[°R, j-e< T i - T e> + °R, ' T e> + G R,k-e< T k ' T e>] + <<C V e> 


( 5 ) 


The change in temperature of fuel-pin node e with respect to time is obtained by re- 
arranging equation (5); that is, 


AT, 


A " ' ) e V e C p,e 


7 — [ G C,f-e< T f T e> + G C,d-e< T d- T e> + G R,j-e( T f - T e> 


G R, l-e^t - T e> + G R,k-e< T k ‘ T e> + Oi" V e] 


( 6 ) 


Now consider a typical node in the honeycomb tube, such as node k in figure 8. The 
thermal resistance network for the honeycomb tube nodal point k is shown in figure 10. 
Applying equation (4) to nodal point k (which has no internal heat source) we have 


AT, 


p k v k c p,kT7 ' G C, ' T k> + G C,j-k< T i - T k> + g r, ;-k< T ; - T k> 


+ G R,f-k (T f ” T k ) + G R,e-k <T e " T k* + G R,d-k* T d ' T k* + G R,j-k* T / ‘ T k l 


G R,i-k< T i 4 - T k> + G R,g-k( T g- T k) 


( 7 ) 


14 






Figure 10. -Thermal resistance network for honeycomb tube nodal point k. 


Similarly, the change in temperature of node k with respect to time is 


AT, 




■7 — [ G c, l-k (T l ■ T k> + G C, j-k <T j ■ T k ) + g R, l-k (T l ' T k ) 


+ G R,f-k< T f- T k> + G R,e-k< T e- T k> + G R,d-k< T d- T k> + G R,j-k< T , 4 ‘ T k> 


G R,i-k( T i 4 - T k> + G R,g-k( T g- T Q 


( 8 ) 


Relations similar to equations (6) and (8) describe the thermal response of each node 
in the 30° sector model of the core. 

The CINDA-3G program computes the temperature change with respect to time of 
each node in the thermal network. As stated in the section Analytical Model of Core, 
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we used the CINDA-3G forward difference subroutine called CNFRWD to determine the 
time rate of change in temperature of each node in the network. The size of the time 
step (i.e. , the value of A 0) used in the simultaneous solution of the thermal network 
equations is defined by a stability criteria formulation in the CNFRWD subroutine. For 
the time period A 0, the temperature change with respect to time of any node is 


AT 

A0 



(9) 


And the temperature of any node at time 0 2 is related to the node temperature at the 
preceding time, 0^, by the expression 

T a = T fl + AT (10) 

e 2 e x 


The basic assumption used in this loss- of- coolant analysis is that all heat inter- 
change between surfaces in the core is by thermal radiation. Because thermal radiation 
is the governing mode of heat transfer, the surface emissivity has a significant effect on 
the thermal response of the core. 

The published values of emissivity for the materials in the core show wide vari- 
ances. Because of these reported variations, we considered three different values of 
emissivity for the loss- of- coolant analysis, namely, 0.20, 0.40, and 0.80. 

All radiating surfaces in the core were assumed to have the same emissivity. And 
a separate analysis was made for each of these three values of emissivity. 


RESULTS AND DISCUSSION 

As we have stated, two different schemes were considered for removing decay heat 
from the core following a total loss of reactor coolant. For each of these schemes, we 
determined the temperature response of all nodes in the 30° sector model of the core. 
But a description of the temperature- time history of all nodes is beyond the scope and 
ambitions of this report. Therefore, in what follows, we shall confine our discussion 
to the temperature response of the fuel pins which are located in the main body of the 
core. 

First, we will describe results for the scheme in which the 4 n shield which sur- 
rounds the core is considered as a single, constant- temperature sink for absorption of 
heat transferred from the core. Then, we will describe results for the second scheme, 
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in which a redundant coolant channel at the center of the core is used as a thermal sink 
in conjunction with the surrounding An radiation shield. 


Scheme I - Decay Heat Absorption by Only the Surrounding 4tt Radiation Shield 

For this scheme of decay heat rejection, the direction of heat transfer is radially 
from the interior of the core to the surrounding An radiation shield. The radiation shield 
temperature is assumed to remain constant at 700 K (1260° R) following the loss of 
coolant. 

Temperature response of core centerline fuel pin . - The centerline fuel pin (fuel 
pin 1 of fig. 5) has the largest power factor and also the largest overall resistance to 
radiation heat transfer. Thus, with a complete loss of reactor coolant and with radia- 
tion as the controlling mode of heat transfer, the centerline fuel pin reaches the highest 
temperature of any fuel pin in the core. 

Figure 11 shows the temperature response of the centerline fuel pin for surface 
emissivities of 0.20, 0.40, and 0. 80. From this figure we see that with an emissivity 
of 0.20 the centerline fuel pin reaches the melting point of UN at approximately 1900 sec- 
onds after the loss of coolant. In figure 11, the temperature response curve for the 
0. 20 emissivity is extrapolated (shown dashed) beyond the fuel melting point. It should 
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5700 


.E 5300 1— 
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4500 
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3300 — 


E 2900 U 


25001— 



Figure 11. - Temperature of centerline fuel pin in fast -spectrum, space-power reactor against time after loss 
of coolant, for surface emissivities of 0.20, 0.40, and 0.80. 
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be recognized, however, that the dashed portion of this curve is not realistic. At the 
melting point, the analytical model for heat transfer breaks down due to the latent heat 
effects accompanying a phase change and geometry changes resulting from fuel melting. 
Hence, the dashed portion of the 0. 20 emissivity curve represents only a ficticious tem- 
perature response that would exist if melting did not take place. 

The significance of the data in figure 11 is that with loss of coolant from the core, 
fuel melting occurs only when the surface emissivity is somewhat less than 0. 40. 

Table IV lists the maximum temperature of the centerline fuel pin and the times required 
to reach these maximum temperatures for the three emissivity values. Thus, from 
figure 11 and/or table IV, we conclude that for emissivity values greater than about 
0. 40, thermal radiation to only the surrounding 4ff shield is sufficient to prevent the 
fuel from reaching its melting point. 


TABLE IV. - TEMPERATURE RESPONSE DATA OF CORE 
CENTERLINE FUEL PIN 


Surface 

emissivity 

Maximum temperature 
of centerline fuel pin 

Time to reach melting point 
or maximum temperature, 
sec 

K 

°R 

0.20 

(a) 

(a) 

b 1900 

.40 

2820 

5080 

4000 

.80 

2355 

4240 

2200 


a Greater than fuel (UN) melting temperature, which is approxi- 
2923 K (5260° R). 
xo reach melting point of UN. 


Temperature response of fuel pins in main body of core . - In this section we pre- 
sent a brief description of the temperature response of the individual fuel pins in the 
main body of the core for emissivity values of 0.20, 0. 40, and 0. 80. The numbering 
scheme which defines the fuel pins and nodal points in the analytical model of the core 
is given in figure 5. In the discussion which follows we refer to this numbering scheme 
(fig. 5) to identify the fuel pins and nodal points. 

Emissivity of 0.20: We have shown in the preceding section (fig. 11) that with an 
emissivity of 0. 20, the core centerline fuel pin reaches the melting point of UN at about 
1900 seconds after the loss- of- coolant accident. As we have stated, with fuel melting, 
the analytical model of the core is no longer applicable. Figure 12 shows the tempera- 
ture response of the fuel pins in the main body of the core for the 0.20 emissivity. The 
temperatures in figure 12 are fuel- pin average values. That is, they are the average 
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Fuel -pin average temperature, 







of the nodal point temperatures in each of the respective fuel pins. As a point of inter- 
est, with an emissivity of 0. 20, the maximum temperature difference between modes 
associated with any single fuel pin is small - less than 20 K (36° R). However, the 
circumferential temperature gradient in some of the T-lll honeycomb tubes is consider- 
ably larger. At some locations in the main body of the core, the extreme temperature 
difference between nodes in a single T-lll honeycomb tube is as high as 220 K (396° R). 
For the 0. 20 emissivity, all calculated temperature data beyond 1900 seconds represent 
only the response that would occur if there were no melting in the core. For this rea- 
son, all data curves in figure 12 are shown dashed for times greater than 1900 seconds. 
The spacing between curves in figure 12(a) is a measure of the radial temperature gra- 
dient in the core. From this data, we note that the radial temperature gradient between 
pins is relatively small near the center of the core and increases significantly with 
radial distance from the core center. The radial temperature gradient depicted in fig- 
ure 12(a) is a characteristic of this thermal network in which the resistance to radiation 
heat transfer between the fuel pins increases with decreasing fuel- pin temperatures. 

Emissivity of 0. 40: Figure 13 shows the calculated temperature response of the 
fuel pins in the main body of the core based on an emissivity of 0. 40. From figure 13, 
we note that the temperatures of all fuel pins in the main body of the core are below the 
melting point of UN. However, some of the fuel pins (e.g. , fuel pins 1, 2, 3, and 4) 
reach temperatures which are within about 200 K (360° R) of the UN melting tempera- 
ture. Hence, with an emissivity of 0. 40, the margin between the maximum fuel- pin 
temperature and the melting point of UN is relatively slim, particularly for those pins 
in the center region of the core. 

Emissivity of 0. 80: With an emissivity of 0. 80, the temperatures of all fuel pins 
in the main body of the core are considerably below the UN melting temperature. Fig- 
ure 14 shows the calculated temperature response of these individual fuel pins based on 
an emissivity of 0. 80. In figure 14, the difference between the highest fuel-pin temper- 
ature (that of fuel pin 1) and the melting temperature of UN is approximately 570 K 
(1026° R). Thus, with an emissivity equal to 0.80, there is a significant margin between 
the maximum fuel- pin temperatures and the melting temperature of UN. 

In the foregoing section, we have shown the effect of emissivity on the temperature 
response of the individual fuel pins when only the 4zr shield is considered as a heat sink. 
With this scheme of removing decay heat, fuel melting takes place only when the emis- 
sivity is somewhat less than 0. 40. That is, for surface emissivities greater than about 
0. 40, thermal radiation to the surrounding 47 t shield is sufficient to prevent the fuel 
from reaching its melting temperature. 

Heat interchange between reactor and shield . - Figure 15 shows the total rate of 
heat transfer from the reactor to the surrounding radiation shield as a function of time 
after loss of coolant. For each emissivity, the rate of heat transfer decreases imme- 
diately following the loss of coolant, reaches a minimum, and then increases with time 
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Figure 15. - Total rate of heat transfer from core to shield against time after loss 
of coolant for surface emissivities of 0. 20, 0. 40, and 0. 80. 


and approaches a relative maximum. Also shown in figure 15 is a dashed curve which 
represents the total rate of decay heat generation for this reactor against time after the 
loss of coolant. The difference between the total decay heat generation rate and the total 
rate of heat transfer from the core is the rate of thermal energy storage in the core. 

A significant point of figure 15 is that the thermal capacitance of the core (and hence 
the average core temperature) continues to increase until the total decay power is equal 
to the total rate of heat transfer from the core, that is, until the respective curves 
intersect. Figure 15 shows that with an emissivity of 0. 80, the curves intersect at 
about 4200 seconds. But with emissivities of 0.40 and 0.20, intersections occur at 
times greater than 7000 seconds. Thus, with the 4 n radiation shield as the only sink for 
heat removal from the core, a large fraction of the decay power goes toward increasing 
the thermal capacitance of the core. 
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Scheme II - Decay Heat Absorption by Both a Redundant 
Coolant Channel and the Surrounding 47r Shield 


In this scheme, decay heat generated in the core is transferred to both the redun- 
dant coolant channel at the center of the core and the radiation shield which surrounds 
the core. The redundant coolant channel has a outside diameter of 2. 16 centimeters 
(0. 85 in. ), which is the same as that of the honeycomb tubes. For this scheme, we 
used a constant sink temperature of 1222 K (2200° R) for the core centerline coolant 
passage and a constant sink temperature of 700 K (1260° R) for the radiation shield. 

Temperature response of fuel pins in main body of core . - Figure 16 shows the 
temperature response of the fuel pins in the main body of the core for this scheme of 
decay heat removal. The data shown in this figure were calculated based on an emis- 
sivity value of 0. 20. As indicated in figure 16, the temperatures of all fuel pins in the 
main body of the core are below the melting point of UN. The highest fuel temperature 
(that of fuel pin 5) is approximately 2840 K (5100° R), or about 90 K (160° R) less than 
the melting temperature of UN. 

A comparison of the data in figures 16 and 12 shows that the addition of a redundant 
cooling channel at the center of the core results in a substantial decrease in the fuel- pin 
maximum temperatures. Even with an emissivity value as low as 0. 20, the tempera- 
tures of all fuel pins remain below the melting temperature of UN. 

Heat-transfer rates to centerline cooling passage and surrounding 4n shield . - Fig- 
ure 17 shows the rates of heat transfer to both the centerline cooling passage and the 
surrounding 4 tt shield against time after loss of coolant. The total rate of heat removal 
from the core is the sum of the individual rates represented by the solid curves in this 
figure. 

The dashed curve in figure 17 represents the total rate of decay heat generation in 
the reactor. From this figure, we see that the rate of heat removal by the centerline 
cooling channel is small in relation to the total decay power. Nevertheless, it is suffi- 
cient to prevent the temperatures of the fuel pins from reaching the melting tempera- 
ture of UN, even with an emissivity value of 0.20. 
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Figure 17. - Total rates of heat input to core centerline cooling passage and surrounding 
4 tt shield against time after loss of reactor coolant. Emissivity, 0. 20. 


SUMMARY OF RESULTS 

A two-dimensional, transient, heat- transfer analysis was made to determine the 
temperature response in the core of a conceptual space-power nuclear reactor following 
a total loss of reactor coolant. For this analysis, we assumed that (1) the reactor was 
operative continuously at design power (2. 17 MW) for 1 year prior to the loss of coolant 
and (2) the reactor was shut down immediately following the loss of coolant. With loss 
of coolant from the reactor, the controlling mode of heat transfer is thermal radiation. 

In one of the schemes considered for removing decay heat from the core, we 
assumed that the 4?r shield which surrounds the core acts as a constant- temperature 
sink (temperature, 700 K (1260° R)) for absorption of thermal radiation from the core. 
Results based on this scheme of decay heat removal show that the surface emissivity 
has a significant effect on the temperature response of the core. With an emissivity of 
0. 20, the core centerline fuel pin reached the melting temperature of UN (2923 K, or 
5260° R) at about 1900 seconds after the loss- of- coolant accident. And with emissivity 
values of 0. 40 and 0. 80, the maximum fuel- pin temperatures were 2820 K (5080° R) and 
2355 K (4240° R), respectively. Thus for emissivity values greater than about 0. 40, 
thermal radiation to only the surrounding 4 jt shield is sufficient to prevent the fuel from 
reaching its melting point. 

In another scheme for removing decay heat, the centerline fuel pin was replaced by 
a redundant coolant channel. With this arrangement, decay heat generated in the core 
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is transferred to both the redundant coolant channel at the center of the core and the 4? r 
radiation shield which surrounds the core. A constant temperature of 1222 K (2200° R) 
was assumed for the redundant coolant channel, and the An shield temperature was 
assumed constant at 700 K (1260° R). 

Based on a surface emissivity of 0.20 for all core materials, the calculated maxi- 
mum fuel temperature for this scheme of heat removal was 2840 K (5100° R), or about 
90 K (160° R) less than the melting temperature of UN. Hence, it appears that a single 
coolant channel at the center of the core is sufficient to prevent fuel melting, even with 
a surface emissivity value as low as 0.20. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 8, 1971, 

112-27. 


32 


APPENDIX - SYMBOLS 


o p 

A cross-sectional area for conduction, m , ft 

2 2 

A g surface area, m ; ft 

Cp specific heat, kW-sec/(kg)(K); Btu/(lbm)(°R) 

F configuration factor for radiation heat transfer 

Gq conductor conductance, kW/K; Btu/(sec)(°R) 

°R radiation conductance, kW/K 4 ; Btu/(sec)(°R) 4 

k thermal conductivity, kW/(m)(K); Btu/(sec)(ft)(°R) 

l characteristic length between conductor nodes, m; ft 

Q O 

q"’ volumetric heat generation rate, kW/m ; Btu/(sec)(ft ) 

T temperature, K; °R 

V volume, m 2 ; ft 2 

e emissivity 

6 time, sec 

p density, kg/m^; lbm/ft 2 

ct Stefan- Boltzmann constant, 5.71x10"^ kW/(m 2 )(K 4 ); 0. 476x10” Btu/(sec) 

(ft 2 )(°R) 4 

Subscripts: 

nodal designations (fig. 5 and/or fig. 8) 

general nodal point designation 
general nodal point designation 
time at beginning of time period - 0j) 
time at end of time period (fig - 0j) 


g, i, J, > 

k ,1 J 

a 

0 

e l 

e 2 
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